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Summary 

The  dielectric  properties  of  polar  solvents  are  considered  within  the  context  of  generic  mean 
spherical  approximations.  Given  the  molecular  diameter  in  a  spherical  representation,  the  dipole 
moment,  and  molecular  polarizability,  the  static  permittivity  of  such  a  solvent  is  estimated  on  the 
basis  of  the  Wertheim  polarization  parameter,  and  a  stickiness  parameter  which  accounts  for  other 
electrostatic  features  of  the  bulk  system  and  non-sphericity  of  the  actual  solvent  molecules.  The 
stickiness  parameter  is  related  to  the  magnitude  of  dipole-dipole  interactions  for  a  series  of  solvents 
with  the  same  polar  group.  This  model  is  also  used  to  estimate  the  Kirkwood  structure  factor.  The 
values  obtained  are  in  excellent  agreement  with  structural  features  of  these  systems  obtained  from 
experiment. 
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Introduction 

In  a  previous  paper^  (referred  to  as  I),  we  discussed  the  mean  spherical  approximation 
(MSA)  with  sticky  association  in  a  model  for  the  dielectric  properties  of  polar  solvents.  In  this 
model,  the  molecules  are  represented  as  hard  spheres  with  a  permanent  dipole  moment,  a 
polarizability,  and  an  anisotropic  stickiness  parameter  which  leads  to  structure  formation  or  break 
up  in  the  liquid  system.  The  idea  of  representing  the  complex  molecules  of  which  polar  solvents 
are  composed  as  spheres  is  useful  for  a  large  number  of  applications.  The  analytical  solution  of  the 
MSA  for  polar^  and  polar-polarizable  hard  spheres3.4  has  opened  up  new  perspectives  for  the 
simple  representation  of  these  rather  complex  systems. 

Water  in  its  liquid  state  has,  presumably  the  local  structure  of  ice  l5>6,  in  which  a  central 
water  molecule  is  surrounded  by  four  others  in  tetrahedral  positions.  Recently,  we  have  shown^ 
that  an  analytical  potential  with  terms  up  to  the  octupole  can  represent  very  well  the  experimental 
atom-atom  pair  correlation  functions.  All  of  these  works  are  based  on  the  invariant  expansion 
formalisms.^  in  which  the  pair  (correlation,  potential)  functions  are  expanded  in  terms  of 
rotationally  invariant  generalized  spherical  harmonics. 

Picking  a  reference  frame  in  which  the  intermolecular  center  to  center  vector  is  the  z-axis,  one 
can  classify  the  harmonics  by  the  periodicity  of  the  functions  for  rotations  around  that 
intermolecular  axis,  very  much  like  the  so-called  Hylleras  molecular  orbitals  for  the  H2+  radical: 
the  lowest  term  is  the  x  =  0  (2  orbital)  coefficient^,  which  corresponds  to  the  cylindrically 
symmetric  representation;  the  next  %  =  1  representation  has  the  symmetry  of  cos  (]),  namely  1  node 
for  a  full  rotation  0  <  <j)  <  2jt  around  the  z  (intermolecular)  axis.  Only  the  x  =  0  and  x  =  1 
representations  are  necessary  to  obtain  an  almost  quantitative  representations  of  the  structural  pair 
correlation  functions  of  water.  10 

The  same  representations  are  required  to  formulate  the  theory  of  the  dielectric  constant  for 
polar-polarizable  molecules.^  Therefore,  the  Blum-Fawcett  model in  which  solvents  are 
represented  by  hard  spheres  with  a  dipolar  polarizability,  and  an  orientational  sticky  parameter  can 
be  thought  as  being  a  mean  field  theory  for  the  full  molecule,  which  admittedly  has  many  more 
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directional  interactions.  The  parameter  to  of  the  sticky  interaction  is  a  sort  of  average  of  the  effects 
of  the  very  directional  hydrogen  bonds  in  water,  which  have  the  effect  of  aligning  (to  >  0)  or 
misaligning  (to  <  0)  the  hydrogen  bonds.  In  a  way  it  is  a  realization  of  the  structure  forming 

structure  breaking  parameter  of  Frank  and  Wen.'^ 

In  I,  the  static  dielectric  permittivity  of  the  polar  liquid  Eg  was  related  to  the  Wertheim 

polarization  parameter  fusing  the  equation^ 

(1+X/)^  =  16  Eg 

This  relationship  is  not  quite  correct  in  the  presence  of  sticky  directional  interactions,  and  a  new 
relationship  is  derived  below.  As  a  result,  the  parameterization  for  the  solvents  derived  earlier^ 
changes  so  that  a  new  set  of  parameters  are  presented  in  this  paper. 

In  the  following  section,  we  present  a  revised  theory  for  the  dielectric  permittivity  of  a  polar 
solvent.  Then,  the  Kirkwood  gK  factor  is  derived  and  its  relation  to  older  methods  of  estimating  it 
is  presented. 

In  the  generic  MSA  (GEMSA)11>12  approximations,  the  MSA  closure  relation  is 

c(r)  =  -Pu(r)  for  r  >  A » a  (2) 

where  c(r)  is  the  direct  correlation  function,  u(r)  the  potential  energy  for  dipole-dipole  interactions, 
a,  the  molecular  diameter.  A,  a  large  distance  which  describes  a  sphere  containing  the  system,  and 
P  =  1/kT.  It  has  been  shown  that  the  GEMSA  approximations  all  share  the  property  that  the 
Kirkwood  gK  factor  is  internally  consistent  in  the  sense  that  the  different  ways  of  computing  it 
yield  the  same  result.  On  the  other  hand,  this  is  not  necessarily  true  for  other  approximations.  An 
excellent  and  comprehensive  discussion  of  the  dielectric  behaviour  and  its  molecular  origins  can  be 
found  in  the  review  by  Stell,  Patey  and  H0ye.  V arious  methods  of  extracting  the  gK  factor  and 
the  infinite  frequency  limits  are  also  considered  in  the  following  section. 

The  Dielectric  Constant  for  the  GEMSA 

Consider  a  system  of  molecules  with  dipolar  symmetry.  The  correlation  and  potential 
functions  can  be  expanded  using  the  invariant  expansion.2.8-9  in  general. 


f(l,2)  =  f000(r,2) +  <I)ii0fii0(ri2)  +  fn2(i-j2) 


(3) 
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where  1  is  the  shorthand  for  (ri,  0i,  hi  the  position  and  orientation  of  molecule  land  2  is 
shorthand  for  (t2,  62,  <t>2)>  those  of  molecule  2.  rn  =  1  ri  -  r2 1  is  the  center  to  center  distance. 
The  invariants  Otto  and  Ott2  are  scalar  products  of  the  dipole  vectors  pi  and  pi- 

Ott0ocpj'-p2 

Ott2  oc  3(pi  •  ri2)  (P2  •  fn)  -  Pi  •  P2 
Using  the  trigonometric  relationship 

Pi  •  P2  =  cos  Yi2  =  [cosBi  COS02  +  sin0i  sin02  cos  ((l)i-(t)2)]  (6) 

the  irreducible  invariants.  Ox ,  used  below  are  obtained: 


Oo  =  Otto+Ott2  a  COS01  COS02 


(7) 


and 


oV  +  ^-1  =  ■  2  Otto  a  sin0i  sin02cos  ((l)i-(t)2)  (8) 

Consider  a  system  of  roughly  spherical  molecules  which  have,  in  addition  to  short  range 
interactions,  long  ranged  electrostatic  interactions  of  the  type^’^ 


u(l,2)  =  -5"  [Pl  P2'  3  (PI  '12)  (P2' '■12)]  £1^4,112  (9) 

'12  ''3 

In  the  notation  used  here,  the  vector  ri2=  {ri2i  ri2}  with  magmtude  ri2  has  a  direction  given  by 
the  unit  vector  {^12}  =  {012,4>12}.  The  GEMS  A  approximations  require  that  the  direct 
correlation  function  satisfies  the  equation 

(10) 


c(l,2)=  -pu(l,2) 


for  ri2  >  A 


Since  the  dipole-dipole  interaction  is  long  ranged,  the  u(l,  2)  interaction  yields  a  virial  coefficient 
that  diverges  as  the  volume  of  the  sample  goes  to  infinity.  This  is  also  true  for  the  integrals  of  both 
the  pair  correlation  function  g(l,  2)  and  the  direct  correlation  function.^.!  l  Therefore,  the  analysis 
of  the  dielectric  properties  of  fluids  must  include  the  proper  treatment  of  the  long  ranged  dipolar 
contributions  in  the  correlation  functions. 


Consider  first  the  simplest  expression  for  the  Kirkwood  factor  gK- 
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£s-l 

£s+2 


'v 

a 

|p2  +  - 

r 

d(2)  G?(2)p(l).p(2)  ' 

[  Ps-^ 

> 

(11) 


Using  the  notation  of  Nienhuis  and  Deutchl^,  is  the  number  density  of  dipoles,  and  Q,  a  solid 

angle.  The  function  G2  (2)  =  g(l,  2)  is  the  pair  correlation  function  for  the  system  in  an  external 

electrostatic  test  field  Eg.  The  pair  correlation  function  has  an  angular  dependence.  One  may  show 
that  a  completely  general  expansion  of  this  function  is 

02(2)  =  g(l,2)  =g000(ri2)  +  g"®(ri2)-^J^  +  (Uz)  (1^) 


The  integral  in  equation  (11)  contains  all  three  terms  of  the  expansion.  The  first  term  is  zero 
because  of  the  angular  integrations.  The  second  term  yields  the  Kirkwood  gK  factor: 


gK 


1  +  ditps 


J  drgii®(r)r2 

0 


(13) 


In  the  older  literamre,  the  integral  containing  the  last  term,  gii2  was  neglected  or  improperly 
estimated.  If  the  volume  of  the  system  were  tmly  infinite,  this  integral  would  vamsh  because  of 
orthogonality  in  the  angular  contributions,  and  diverge  because  of  the  radial  term. 

Using  the  notation  of  previous  work2.3.l  1.12,15^  equation  (1 1)  may  be  rewritten  as 


si  ^  Z  [  1  +  4n  ps  o'  (k  =  0)  +  2(1+  47tps  (k  =  0))] 

es+2  3 


l-47tps  cV  (k=0)  l-47Cps  cV  (k=0)J 


(14) 


where 


4jt  p  pp^ 
9 


y 


(15) 


This  result  is  incorrect  because  it  does  not  include  the  long  ranged  interactions.  When  these  are 
included  following  Wertheim^'^,  the  following  result  is  obtained 
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es-1 

1 

r  1  ,21 

(16) 

8s+2 

"  3 

Lyo-2y  yi+yJ 

where 

2 

Yo  = 

1  -  4k  Ps  C  0  (k  —  0) 

(17) 

and 

2 

Yi  = 

1 

oo 

1! 

O 

(18) 

The  closure  of  the  GEMS  A  is 

12  2 
y  =  3  (Vo  -  vi) 

(19) 

Applying  this  result  to  equation  (14),  the  result  is 

2  2 

Es-i  _  Vo  -  yi  (20) 

es+2  Yo  +  2y  1 
from  which  one  obtains 

ri-47c  ps  Co\k  =  0)]  _  Yo 

es  = - 1  “2  ^  ^ 

[i-47t  psCii(k  =  0)]  yi 


This  is  an  extension  for  our  case  of  Wertheim's  result.2-4 

Now  explicit  relationships  for  the  anisotropic  sticky  hard  dipole  model  can  be  written.  These 


are 


2  ~ii  X^(k+2)^ 

yo  =  l-47ipsCo  (k  =  0)  =  -^9 - 


1  . 


and 


y  1  —  1  "  47tps  C  o  (k  —  0)  — 


16  XHX+2)'^ 


(22) 


(23) 


>.2(i+)t)4  9 

where  to  is  the  longitudinal  sticky  factor  defined  in  our  previous  work.l  Equations  (21)  -  (23) 
immediately  lead  to  the  relation 
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?l2(l+X)4  to^2  (24) 

'==—16—  (‘-xj 

Furthermore,  we  note  that  the  closure  relation  for  the  sticky  MSA  is 


dp  =  3y  =  (l-to/?^)2(l-l/es)  =  y^yi 

We  consider  now  the  polar-polarizable  system  which  was  investigated  in  detail  by 
Wertheim .3  The  added  complication  for  such  a  system  is  that  the  polanzation  force  is  a  long  range 
many  body  force  and  collective  effects  cannot  be  ignored.  In  other  words,  a  simple  pairwise 
additive  potential  is  not  adequate.  Therefore,  a  renormalization  approximation  with  an  effective 
dipole  moment  that  depends  not  only  on  the  polarizability  but  also  on  the  density  and  molecular 
interactions  is  needed.  A  simple  but  powerful  method  for  computing  these  effects  was  described 
by  H0ye  and  Stell.16  it  is  based  on  the  concept  that  a  system  of  polar-polarizable  molecules  is 
equivalent  to  a  mixture  of  molecules  with  permanent  dipole  moments  ma  where 

iHa  =  Pa  Ep  ■  ^ 


where  a  is  the  polarizability  tensor  and  Ep  is  the  local  (fluctuating)  microfield.17  ft  can  be  shown 
that  this  yields  a  random  Gaussian  distribution  of  the  electric  microfield,  and  therefore,  of  the 
instantaneous  dipole  moment.  However,  as  shown  by  Adelman  and  Deutchl^,  a  random  imxture 
of  equal  size  molecules  of  diameter  as  with  fluctuating  dipoles  in  the  MSA  and  the  GEMSAs  is 

equivalent  to  a  fluid  with  a  single  dipole  moment  given  by 

<p>  =  E  Pi  <mi>2/  E  Pi 
The  result  for  the  MSA  approximation  is 


where 


2 

d2  +  47t  a  ps  ot 


1 

8a  A-1  ^ 

Os 


(28) 

(29) 
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and 

^2  ^  4k  PPs  (30) 

This  gives  different  results  from  our  previous  analysis  because  the  factor  (1  -  to  /  is  now 
included.  Finally,  the  Kirkwood  gK  factor  can  be  defined  in  several  different  ways,  all  of  which 
are  equivalent  in  the  GEMSA  type  of  approximations.  Thus,  one  obtains 


gK=  1  +Psh"<>(0) 

or 

(£s-1)(2£s+1) 

gK  =  2 

3£s  dp 


(31) 

(32) 


Experimental  Data  for  Polar  Solvents 

In  I,  experimental  permittivity  data  at  both  static  and  optical  frequencies  were  related  to 
molecular  properties  of  polar  solvents  using  the  MSA.  In  this  regard  the  solvent  molecules  were 
assumed  to  be  approximately  spherical  with  diameters  estimated  from  gas  phase  solubilities.  19  The 
molecular  diameter  as  together  with  other  molecular  properties  are  given  in  Table  1  for  20  polar 

solvents,  both  protic  and  aprotic. 

The  molecular  polarizability  may  be  calculated  from  the  optical  permittivity  using  the  present 
model  as  described  in  I.  The  Wertheim  equation  is  then  written  as 

eop  =  (^.p+D-'/W 

where  e<,p  is  the  square  of  the  refractive  index,  and  Xpp,  the  value  of  the  polarization  parameter  at 
optical  frequencies.  One  may  then  estimate  the  optical  value  of  the  parameter  dp  using  the 
relationship 

2 

,2  _  )lop  (X,op+2)2  .  J_N  (34) 

=  9  I  e„pj 

Alternatively,  dop  can  be  calculated  from  equation  (28)  written  in  a  form  applicable  at  optical 


frequencies: 
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2 

dop 


47C  a  ps  oc 


where 


1 


^  ^  /^Xpp-l  ^ 

<^5^  V^op+2 ) 

The  polarizability  a  is  found  by  solving  equation  (35) 
equation  (34). 


(35) 

(36) 


using  the  value  of  dop  obtained  from 


Values  of  a  found  by  this  route  are  also  recorded  in  Table  1.  These  estimates  are  slighdy 
higher  than  those  obtained  from  the  classical  continuum  model.20  Equations  (33)  -  (36)  represent 
an  approximate  solution  of  the  high  frequency  problem.  However,  since  there  is  good  agreement 
between  values  of  a  obtained  experimentally,  from  the  continuum  model  and  from  the  MSA 
model,  we  believe  that  the  present  treatment  deals  adequately  with  the  contribution  to  the 
permittivity  from  molecular  polarizability  at  static  frequencies. 

Estimation  of  the  MSA  parameters  at  low  frequencies  involves  solving  the  following  three 

non-separable  equations  for  X  and  to: 


?i2(l+?.)4 

16 


(37) 


2  _  X2(U2)2  r.  t^Y 

'p  "  9  Li  x> 


and 


2 

dp 


a2  d2  +  47C  a  Ps  a 


(38) 


(39) 


Equation  (38)  is  obtained  by  subtracting  equation  (23)  from  equation  (22).  In  order  to  obtain  a 
solution,  to  was  initially  assumed  to  be  zero.  An  initial  estimate  of  X  was  then  obtained  from 
equation  (37).  Using  the  solvent's  molecular  density,  dipole  moment,  and  polarizability,  the 
parameter  dp  was  estimated  on  the  basis  of  equation  (39).  Finally,  the  stickiness  parameter  to  was 
estimated  using  equation  (38).  This  process  was  repeated  using  the  new  estimate  of  to  in  equation 
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(37).  Eventually  by  iteration,  values  of  X  and  to  corresponding  to  an  exact  solution  of  these 
equations  were  obtained. 

Values  of  the  MSA  parameters  determined  for  22  polar  solvents  are  summarized  in  Table  I. 
the  numerical  values  of  X  and  to  are  different  than  those  obtained  earlierl  but  follow  the  same 
qualitative  trends.  In  the  case  of  aprotic  solvents,  to  is  negative.  This  indicates  that  the  polarization 
parameter  X  is  overestimated  when  the  simpler  model  based  on  non-sticky  interactions  is  used.  In 
addition,  these  interactions  result  in  the  break  up  of  structure  due  to  head-to-tail  interactions. 

Another  trend  apparent  in  the  data  in  Table  I  shows  that  parameter  to  increases  in  magnitude 
as  the  length  of  the  alkyl  group  increases  for  solvents  with  the  same  functional  group.  Thus,  in  the 
nitriles,  to  increases  from  -2.53  for  acetonitrile  to  -1.93  for  butyronitrile.  Similarly,  for  the 
alcohols,  to  increases  from  0.81  for  methanol  to  1.47  for  1-butanol.  The  positive  value  of  to  in 
these  solvents  indicates  that  head-to-tail  interactions  provide  a  major  contribution  to  the  polarity  of 
the  solvent.  A  similar  comment  can  be  made  about  N-methylformamide,  the  only  other  solvent 
with  a  positive  value  of  to- 

As  noted  earlierl,  the  stickiness  parameter  is  linearly  related  to  the  strength  of  dipole-dipole 
interactions  within  a  given  series  of  solvents.  This  is  illustrated  in  Figure  1  for  the  water-alcohol 
and  alkyl  nitrile  series.  As  the  length  of  the  alkyl  chain  increases,  the  parameter  to  becomes  more 
positive  as  noted  above.  This  trend  accompanies  the  displacement  of  the  dipolar  functional  group 
further  from  the  geometrical  center  of  the  molecule.  At  the  same  time,  the  tendency  for  structure 
forming  interactions  in  head-to-tail  configurations  increases. 

Values  of  the  Kirkwood  structure  factor  gK  calculated  according  to  the  MSA  (equation  (30)) 
are  also  given  in  Table  I.  These  parameters  vary  from  1.09  for  butyronitrile  to  4.52  for  N- 
methylformamide.  The  latter  solvent  is  known  to  be  highly  structured  and  to  form  long  chains  of 
solvent  molecules  linked  through  hydrogen  bonding.^l  xbis  results  in  the  very  high  static 
permittivity  and  value  of  gK-  In  the  case  of  water,  the  value  of  gK  is  2.79,  that  is,  very  close  to  the 
value  estimated  by  Kirkwood^^’^^.on  the  basis  of  the  well  known  tetrahedral  structure  for 
surrounding  hydrogen  bonded  water  dipoles.  The  value  of  gK  for  the  alcohols  is  higher. 
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increasing  from  2.99  for  methanol  to  3.26  for  n-butanol.  This  suggests  that  the  tendency  to  form 
hydrogen  bonded  chains  in  these  solvents  increases  with  the  length  of  the  alkyl  group. 

Values  of  gK  for  the  aprotic  solvents  lie  between  1  and  2.  Acetonitrile  is  an  example  of  a 
solvent  with  very  little  structure  in  the  bulk.24  Its  value  of  gK  is  close  to  unity  as  expected  for  n- 
butanol.  This  suggests  that  the  tendency  to  form  hydrogen  bonded  chains  in  these  solvents 
increases  with  the  length  of  the  alkyl  group. 

Values  of  gK  for  the  aprotic  solvents  lie  between  1  and  2.  Acetonitrile  is  an  example  of  a 
solvent  with  very  little  structure  in  the  bulk.23  Its  value  of  gK  is  close  to  unity  as  expected.  On  the 
other  hand,  dimethylsulfoxide  which  is  known  to  be  strongly  associated  as  dimers  in  the  bulk25 
has  a  value  of  gK  equal  to  1.67.  In  the  case  of  propylene  carbonate,  gK  is  equal  to  1.86  suggesting 
some  degree  of  association  in  the  pure  liquid.  In  fact,  dielectric  relaxation  measurements  for  this 
system26-28  show  the  existence  of  two  relaxation  processes  which  could  be  interpreted  as  evidence 
of  molecular  association.  More  than  one  relaxation  process  is  observed  for  several  aprotic 
solvents.29  When  the  dipole  moment  is  high,  dimer  formation  may  take  place  in  these  systems. 

The  values  of  gK  reported  here  are  very  different  from  those  estimated  by  Marcus.39  xhe 
latter  quantities  were  calculated  on  the  basis  of  a  revised  version  of  the  Kirkwood  equation 
proposed  by  Cole.31  Cole's  relationship  is  not  consistent  with  the  dielectric  theory  used  herel2 
and  therefore  is  not  considered  to  give  a  good  representation  of  the  relative  values  of  the  structure 
factors  for  the  solvents  discussed. 

Discussion 

The  revised  formulation  of  the  MSA  presented  here  leads  to  new  values  of  the  parameters  X 
and  to-  However,  it  is  clear  that  the  new  parameters  follow  the  same  trends  as  those  presented 
earlier.^  Thus,  the  conclusions  reached  earlier  regarding  the  stickiness  parameter  to  are 
qualitatively  correct.  The  MSA  model  is  both  analytical  and  relatively  simple,  providing  an 
excellent  molecular  description  of  the  dielectric  properties  of  the  polar  solvents  considered. 

In  applying  any  molecular  model  one  must  have  accurate  values  of  the  molecular  parameters 
involved.  In  the  present  case,  these  are  the  dipole  moment,  polarizability,  and  molecular  diameter. 
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The  diameters  which  are  based  on  gas-phase  solubilities  19  correlate  well  with  the  cube  root  of  the 
molar  volume  and  therefore  are  regarded  as  giving  good  estimates  within  the  context  of  the 
spherical  approximation  used.  In  previous  work  dipole  moments  were  taken  from  compilations 
based  largely  on  liquid  solution  measurements.^^  The  values  in  Table  1  are  mainly  those  from  gas 
phase  measurements^^  and  are  considered  to  be  more  appropriate  to  the  present  correlations.  This 
is  apparent  with  the  revised  dipole  moment  data  for  the  nitriles.  For  these  solvents,  the  values  of  to 
including  that  for  benzonitrile,  as  linear  in  the  energy  parameter  pVd^.  Obviously,  small  errors  in 
p  can  result  in  large  errors  in  p^/d^,  so  that  in  the  previous  correlation,  the  benzonitrile  data  did  not 
fit  with  the  results  for  the  alkyl  nitriles. 

Finally,  the  present  estimates  of  the  structure  factor  gK  appear  to  be  very  reasonable  on  the 
basis  of  what  is  known  about  the  structure  of  the  solvents  considered.  More  experimental  work  in 
this  area  is  definitely  needed.  Infrared  spectroscopy  has  proven  especially  effective  in  the  case  of 
aprotic  solvents  and  provides  a  means  of  understanding  the  estimates  of  gK  obtained  in  these 

solvents.^^’^^ 
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Figure  Legend 

Figure  1.  Plot  of  the  stickiness  parameter,  to  versus  the  dipole-dipole  interaction  energy  factor 
p2/d3  for  the  lower  alcohols  and  water,  and  for  the  nitriles.  The  abbreviations  for  the 
solvents  are  defined  in  Table  1 . 


